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HYDRAULICS OF WELLS* 


Dean F. Peterson, Jr.,} M. ASCE 


INTRODUCTION 


Groundwater is extensively used, not only for irrigation in the Western 
United States and in other arid and semi-arid regions, but all over the world 
for domestic and industrial purposes. The basis for development of this re- 
source is a long history of practical experience plus relatively recent 
developments of theory. 

The study of well hydraulics is a fascinating subject and has occupied the 
attention of many capable investigators. Like all engineering sciences, how- 
ever, it is not an exact science and the mathematical assumptions which often 
have to be made represent only limiting cases of real conditions. For in- 
stance, one interested in theoretical well hydraulics usually finds he must 
postulate at least some of the following conditions, which are not usually fully 
correct: 

1) The natural groundwater regime affecting an aquifer remains 

constant in time. 

2) The permeability is uniform and isotropic throughout the 

formation. 

3) The porosity and thickness of the formation is uniform. 

4) The well penetrates fully and the hydraulic head is the same at 

all points on the well boundary. 

5) The initial piezometric surface is level. 

6) Only laminar flow exists in the region affected by the well. 

7) The region from which the well draws is infinite in extent. 

In spite of the necessity for introducing some of these unrealities, excellent 
progress has been achieved by students of theoretical well hydraulics. Useful 
tools for taking into account specific deviations of natural conditions from the 
idealized conditions assumed have also been devised in many instances. Well 
hydraulics, especially in recent years, has contributed much toward a better 
and more economic development of the important groundwater resource. This 
paper will attempt to review the general development of theoretical well hy- 
draulics and to comment on certain aspects of the application of these devel- 
opments to practical engineering problems. 


Flow Formulas 


Equilibrium Case 

Dupuit, 1863, (1) is credited with first combining Darcy’s Law of laminar 
flow through sands and gravels with the statement of continuity in order to 
derive an equation for well discharge. Dupuit assumed complete axial 


* Paper presented to Irrigation Division, ASCE, Salt Lake City, Utah, 
Sept. 10, 1954. 
1. Prof. and Head, Dept. of Civ. Eng., Colorado A & M College, Fort Collins, 


Colo, 
708-1 


symmetry, steady flow through an infinitely extending aquifer, and ignored the 
effects of curvilinear motion. He deduced 


(h2-h,? 
Q= mK, ’ (1) 


where, Q is the discharge, (L377), K is the permeability or, more properly, 
the hydraulic conductivity, (LT-1), and hi, and hg are the elevations of the 
water surface above a level, impermeable bed bounding the lower side of the 
aquifer at radial distances rj and ro. See Fig. 1. All quantities in Eq 1 must 
be expressed in compatible units or the necessary conversion constants 
introduced. 

As stated in Eq 1, Dupuit’s equation implies a well discharging from a 
water-table condition. However, all of the common steady-flow formulas for 
fully penetrating wells stem from the Dupuit concept and may readily be de- 
duced by algebraic manipulation of Eq 1. Considering Eq 1, 


he* -h,? s (hy—-h,)(h, +h,) = (s, -S,)(2 m) , 


where sj and s2 are the drawdowns from the original water table or piezo- 
metric surface and m is the average thickness of the water-transmitting for- 
mation. Direct substitution yields 


2mKm(s,-s,) 


Eq 2 may also be used to describe a well discharging from an artesian aquifer 
of thickness m. 

Common practice of the Geological Survey is to introduce the term trans- 
missibility T into these equations. Transmissibility is the product Km and is 
widely used in characterizing an aquifer. Making this substitution in Eq 2 
yields 


(2) 


Q In 


In Geological Survey practice T is ordinarily expressed in gallons per day per 
foot and Q, in gallons per minute. Length units are in feet. Introducing the 
necessary unit conversion constants into Eq 3 yields 


527.7 Q 
To (4) 
Se 

Equations J through 4 are useful for determining K or T of water-bearing 
formations using pump tests. This procedure is useful in engineering investi- 
gations of dam foundations and land drainage as well as in water supply work. 
These equations are also useful in predicting well yields. Ordinary practice 
is to plot drawdown against log r, then the quantity (s - 89) /logyg r2/rq = 
As is simply the drop in piezometric surface for one log cycle of radii from 
the well center. T is directly proportional to Q and inversely proportional to 
As, that is T = 527.7 Q/As. 
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If Ss} - Sp = D is considered to be the drawdown of the well, rj corresponds 
to the well radius ry; and rg becomes the radius to a point where the func- 
tioning of the well does not disturb the existing piezometric surface, origi- 
nally assumed as level. This latter distance is ordinarily termed the radius 
of influence. Obviously an equilibrium well cannot exist in a finite, level, 
piezometric field. The concept of a radius of influence accordingly implies 
equilibrium between the flow discharged by the well and the flow replenished 
by nature to the region influenced by the well. The radius of influence is 
simply a length parameter describing the size of the influence region in a 
naturally replenishing field. 

Israelsen, Hansen and the writer (8) in 1950 postulated a general functional 
relationship for steady-state wells expressed in the dimensionless form 


Q 


Kr2 =F(D/r,, m/r,, |), (5) 
where I represents the unit natural replenishment. These authors wrote ex- 
pressions for natural seepage flow intercepted by the influence region of a 
well. Eliminating the radius of influence by substituting these expressions 
into the well discharge equations, Eqs 1, 2, and 3, approximate functional re- 
lationships were obtained for fully-penetrating wells in both artesian and un- 
confined aquifers having natural piezometric slopes ip; and for a well in an 
unconfined aquifer replenished vertically at a uniform rate. Details of these 
solutions appear in the reference. For the unconfined case, this procedure 
results in the equation 


Q 


Khe 2.303 ' (6) 


2 K h.* in 
where h, is the original height of the water surface above the impermeable 
stratum, Figure 1; and i, is the natural slope of the water table. For the 


confined case, 


(7) 


For the case of a well in an unconfined aquifer replenished both by vertical 
seepage at the uniform rate qy (volume per unit area per unit of time) and by 
horizontal seepage into the influence zone, a procedure embracing similar 
assumptions to those of Dupuit in deriving Eq 1 results in, approximately, 


7 K (h? - h,?) 
r ns \2 
-3 (7) 


Q= (8) 


In Eq 8, n is the proportion of total discharge Q originating by vertical 
replenishment and r, is the radius of influence. If n = 1, 
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Kh? 1.15) log| | 
Kh,2/\ W/K 2 
Equations 6, 7, and 9 are not in the identical symbolic form of Eq 5. How- 
ever m and D may be defined in terms of hg and hy for Eqs 6 and 9. I is 
represented by i, in Eqs 6 and 7 and by q,/K in Eq 9. 
The relationships expressed by Eqs 6, 7, and 9 are shown graphically by 
Figures 3, 4, 5, and 6. These should be useful in predicting the performance 
of wells for the limiting case when the discharge is in equilibrium with re- 


charge by nature. 
Letting h = he and r = rg, Eq 9 may be rewritten in the form 


= (het het) (10) 


where e is the naperian base. Equation 10 is interesting because it implies 
that if a transformation ry = rie-n/ 2 is introduced, the resulting flow system 
of a vertically replenished aquifer drained by a well of radius ry, will be 
geometrically similar to that of a Dupuit system drained by a well of radius 
ry. Equation 10 is subjected to the same theoretical approximations as Eq 1 
and in addition to the difficulty that the phreatic surface is not actually a 
streamline as for the Dupuit case. In spite of these approximations, Eqs 9 
and 10 are believed to be useful in the design of relief drainage systems using 
pumped wells. For such systems the average vertical replenishment is apt to 
be uniformly distributed over the influence area. 


Non-Equilibrium Case 

In 1935, Theis, 11) utilizing mathematical methods developed in the study 
of heat flow, proposed a means for analyzing the non-steady flow of water into 
a well. Because the steady-flow or equilibrium condition is a limiting case 
in well hydraulics, rather than the rule, this advance has led to much greater 
understanding of well performance and to improved well engineering. It takes 
into account the water drawn from storage in the aquifer,—as the result of 
unwatering in the water table case, or from volume changes associated with 
changes in pressure under pumping in the artesian case. An additional 
aquifer characteristic, the storage coefficient S was introduced. This coef- 
ficient, in common units, equals the volume of water released in one square 
foot of area tributary to the well owing to a unit decrease in piezometric head. 
For the water table case this becomes the specific yield. 

For the case of an infinitely extending aquifer and constant discharge, 


@ 


ant (11) 
r2S/4Tt 


s 


This integral is known as the exponential integral and tables of its values have 
been published. (10) tt has also been called the well function, W(u), where the 
argument u is the lower limit, r2g/4Tt. Figure 7 shows W(u) as a function 
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of u. For constant discharge, u is proportional to r2/t and s is proportional 
to W(u). By plotting observed values of s vs r2/t to the same scale, respec- 
tively, as W(u) and u one obtains a curve similar to that of Fig. 7. This may 
be matched with the curve of Fig. 7 and a specific match point chosen for 
which values of u, W(u), s, and r2/t may be noted. Rewriting Eq 11, 


$= ae W(u) » (12) 


and recalling, 
r’?/t=4Tu/S , (13) 


one may substitute these specific values to solve for S and T for the forma- 
tion. 

The exponential integral may be expressed by an infinite series. For 
small values of u only the first two terms are significant and W(u) = -0.5772 - 
ln u. Substituting in Eq 12 gives 


- 2 
$= 4nmT (in 0.5772) (14) 


Taking into account the value of u from Eq 13 for a well at a particular radial 
distance one can demonstrate from Eq 14 that 


$,-S, = Inge (15) 


or in Geological Survey units, 
Se -S, = log (16) 


In using Eqs 15 and 16, s may be plotted against log t until the points fall on 
a straight line. Data from the straight line portion of the curve are used to 
solve Eq i6 for T, see Fig. 8. By extending the straight line to a value of t, 
tp, when drawdown is zero, one may combine Eqs 13 and 14 to give a solution 
for S, 


S = 2.25 Tto/r® . (17) 


The work of Theis is based on a boundary condition of constant discharge. In 
1952 C. E. Jacob and S. W. Lohman(5) published a solution for the case of 
unsteady flow to a well at constant drawdown as follows. 


us) 


In Eq 18, a= Tt/sr2, and G(q) is an irregular integral, which could be evalu- 
ated only by numerical methods. Values of G(a@) were also published, For 
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Fig. 8. Drawdown ina Well as a Function of Time. 


708-13 


J 


large values of a, G(a) was found to approach 2/W (1/4 @). The solution of 
Jacob and Lohman is particularly applicable to analysis of flowing artesian 
wells. 

The non-equilibrium formula may be applied to wells drawing from water 
table conditions if proper correction is made to allow for the decrease in ef- 
fective thickness ovtag to lowering of the water table. Jacob(4) has shown 
that a correction of s~/2m should be subtracted from the numerical values of 
observed drawdown, s. Failure to apply this correction results in rather con- 
siderable errors in the computed values of T and S. 

Over the past two decades, the non-equilibrium formula has been widely 
used with excellent results. Even where an equilibrium state may be possible, 
at the beginning of pumping the condition is always one of non-equilibrium. 

As time increases and the cone of influence becomes larger, the condition of 
the region near the well more closely resembles an equilibrium state because 
proportionately less of the flow through that part of the region comes from 
storage. Whether some sort of equilibrium may eventually occur or not de- 
pends on the natural hydrology of the aquifer. 

Operation of a well depletes underground storage. Actually this increases 
the capacity of the underground reservoir to utilize available recharge. From 
an ideal point of view, the management of a system of wells should be such as 
to deplete underground storage annually in an amount equal to available re- 
charge. Continual pumping in excess of this amount would lead to eventual 
depletion; pumping less than this amount would mean some of the available 
recharge waters would be lost owing to lack of storage capacity. 


Method of Images 


By using the common mathematical device of assuming a set of images of 
the real condition, well hydraulics may take into account the effect of finite 
geologic boundaries, deviating from the ideal assumed in deriving the various 
formulas. This may best be illustrated by an example, Consider an aquifer 
bounded by an impermeable boundary or aquiclude as shown by Figure 9. 
Assume the aquiclude to be approximately represented by the vertical plane 
A -B. No flow occurs through this plane. A set of images which makes this 
a plane of symmetry will provide a zero hydraulic gradient at the required 
boundary. The real drawdown will be the sum of the drawdowns for the real 
and the image wells in an infinite aquifer. in the event of two parallel bound- 
aries, two planes of symmetry are necessary and this may be achieved only 
by postulating an infinite set of images. The actual correction need be made 
only for the nearest few of the images, the number depending on the degree 
of accuracy desired. 

In the case involving a source of recharge at constant head, a symmetrical, 
but negative, set of images will provide the necessary boundary condition. 


Hydraulic Effects of the Well 


The foregoing discussion has treated the flow of water through the aquifer 
under an energy gradient created by a well. The water must also be trans- 
ferred through the screen and casing or pump column to the point of dis- 
charge. Under some circumstances the energy expended in moving the water 
through the well structure may exceed that used in moving it through the 
aquifer. Better understanding of hydraulic principles involved in this latter 
mechanism should lead to improved well engineering. Some of the 


708-14 


sebpus jo poyjay “6 


apnjoinby 


p 

| 7 

708-15 


considerations which have come to the writer’s attention are outlined in the 
following discussion. 


Specific Capacity 

Engineers have designated the ratio of discharge to drawdown specific 
capacity. If the hydraulic head losses through the screen and casing were 
zero and the time effect of storage depletion were ignored, the discharge of 
an artesian well could be expected to be directly proportional to drawdown. 
This would lead to a constant value of specific capacity corresponding to all 
values of discharge of an artesian well—a condition usually assumed. For 
wells in unconfined aquifers, an increase in drawdown at the same time de- 
creases the effective thickness of the aquifer. Thus, even discounting energy 
losses at the well, the specific capacity would decrease with discharge for 
the water-table case. 

The hydraulic losses through the well cause further eager | 5 of the 
relationship of discharge to drawdown. As mentioned by Jacob, ‘ flow 
through the screen and casing usually occurs in the turbulent regime and re- 
sulting head losses are thus proportional to Q?. Aquifer losses, under con- 
ditions of laminar flow, should be proportional to Q for an artesian well. Thus 
one may write, 


s=BQ + CQ, (19) 


where s is the total drawdown in the well and B and C are constants. 
Equation 19 can be evaluated approximately by test pumping at two different 
discharges, Q; and Qo, and measuring the respective values, sj and s2. Sub- 


stitution successively into Eq 19 provides simultaneous equations in B and C. 
The difficulty with this procedure is that it is based on the assumption of 
steady flow and does not take into account the effect of time depletion of 
storage. In his excellent paper Jacob(3) has outlined a multiple-step draw- 
down test and procedure which provides a solution for the constants B and C 
and, at the same time, furnishes information leading to determination of the 
effective radius of a well. 


Effective Radius 

The effective radius of a well, as it is used in the formulas of flow, may 
not be the same as the radius of the screen or hole, especially for wells in 
uncemented sediments. Development of the well or use of gravel envelopes 
increases the permeability of the formation immediately surrounding the 
casing. This effect is the same as increasing the radius. The effective 
radius is defined by Jacob as “that distance, measured radially from the axis 
of the well, at which the theoretical drawdown based on the logarithmic head 
distribution equals the actual drawdown just outside the screen.” In the ref- 
erence previously cited, Jacob gives a procedure for determining this quan- 
tity using the results of field tests. 


Seepage Face 

For wells in unconfined aquifers, the Dupuit formula implies that the water 
surface outside of the well intersects the water level in the well at the casing. 
Even for a hydrodynamically perfect well this cannot be the case—a finite 
seepage face always exists between the water level in the well and the water 
surface outside. Hansen,(2) in 1949, proposed a relationship of the form 


K fy 
for a well in an unconfined aquifer bounded by a level impermeable stratum 
at the bottom. In Eq 20, hy is the elevation of the water level in the well and 
hg the corresponding elevation outside, both measured from the impermeable 
stratum. Using experimental data available at that time, Hansen was able to 
propose curves graphing the functional relationship of Eq 20. Additional in- 
formation has since been made available by relaxation investigations of 
Yang(12) and analogy investigations of Zee(13) which make possible presenta- 
tion of Figure 10. Figure 10 also shows values of hj15/rw corresponding to 
the values of the other parameters where hj15 is the height of the water level 
at a distance of 115r,. A field test of a well at Mosca, Colorado by Professor 
W. E. Code and the writer in 1954 correlated very well with the curves of this 
figure although the actual value of hg could not be directly measured. The 
results of this test appear as point A on Figure 10. 

The approach velocity in the neighborhood of an equipotential well is not 
constant, Figure 11. At the top of the seepage surface, the approach velocity 
is K. Theoretically this increases with decreasing elevation and is propor- 
tional to the cosecant of the angle of inclination of the streamlines intersecting 
the seepage face. At the singular point, at the elevation of the water level in 
the well, this velocity theoretically becomes infinite but, practically, turbulent 
flow probably occurs. Below the singular point the velocity rapidly decreases 
to an almost constant value. Figure 11 also shows the relationship between 
Dupuit’s curve and the true water surface, 


Well Screens 

Frequently a large part of the energy imparted through a well is expended 
in transferring the water through the screen and pump. For this reason, at- 
tention should be given to the hydraulic performance of the well structure. 
While considerable progress has been made, collection of data, especially in 
the field, is difficult because of rapidly changing flow conditions near the 
well. Laboratory experiments which similate field conditions are expensive 
and arduous. Nevertheless more attention should be given to this important 
aspect of the problem of well hydraulics. 

An important part of the well structure is the screen. Screens of some 
kind are always required except in cemented sediments. They may range 
from rough, haphazard, perforations in a steel casing to highly engineered 
and carefully manufactured screens of specially selected material. The func- 
tion of a screen is to exclude the natural sediments but to allow the greatest 
possible flow of water into the well. The factor of longevity influences the 
choice of screen. 

The hydraulic performance of well screens was ably treated by Petersen, 
Rohwer, and Albertson. (9) Water enters the interior of a screen in the form 
of radial jets at relatively high velocities. The energy of these jets is dis- 
sipated and the flow accelerated in the axial direction. From a theoretical 
consideration of the mechanics involved, these invesitgators deduced 


Ah =2 cosh(CL/D +1 
V¥2g cosh(CL/D- 4 (21) 


where Ah is the hydraulic head loss involved in the screen, V is the final 
average velocity along the screen axis (Q/A where Q is the well discharge 
and A the cross-sectional area of the screen), L is the axial length of the 
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screen and D is the screen diameter. C was defined as the screen coefficient, 


. (22) 


In Eq 22 C,, is the orifice coefficient of discharge applying to the screen 
opening and Ap is the fractional ratio of screen opening to total screen sur- 
face. Equation 21 was confirmed by large scale laboratory experiments using 
specially-made screens with machined openings arid not surrounded by gravel, 
Figure 12, Some coefficients were also determined for certain commercial 
type screens, 

The loss coefficient, 4h in Eq 22, approaches 2 for values of CL/D 

V*/2¢ 

exceeding about 6. Apparently, for a particular screen of a certain diameter, 
increasing L beyond this amount would not result in appreciable decrease in 
screen losses for a particular discharge. The effect of gravel surrounding 
the screen was to reduce the screen coefficient but not appreciably for larger 
sizes of gravel. Where the gravel size was sufficiently small that it effec- 
tively reduced the percentage of open area of screen, the screen coefficient 
was also appreciably reduced. This important consideration needs to be 
taken into account in transferring the implications of Eq 21 to the field. One 
should note that Eq 21 does not imply that increasing the length of screen 
beyond a value corresponding to CL/D = 6 would fail to result in increased 
discharge, but only that the ratio of screen losses to velocity head along the 
screen axis, for a particular discharge, will not be improved by lengthening 


the screen. 
The action of a well screen is such that the greatest flow occurs near the 


discharge point of the screen and flow through the openings decreases with 
distance from this point. In an actual installation, the discharge point would 
correspond to the intake for the pump bowls. In solving the conventional 
equations for flow through aquifers, the well boundary is assumed as a sur- 
face of constant head equal to the water level in the well. The effect of 
screens is to upset this, the head will be least opposite the pump inlet and 
will increase in both directions from this point. 


Gravel Envelopes 


A well screen is in effect a filter, with an additional characteristic of 
structural rigidity. In many cases, gravel is installed around the screen in 
order to improve the filter action. Such an installation is known as a gravel 
envelope or pack, In other cases a well may be developed in natural materi- 
als by vigorous surging so that the finer sediments are removed and pumped 
out of the well, leaving a natural filter of the coarser materials. 

As for a screen, a desirable gravel envelope excludes the sediments of the 
natural formation and at the same time permits the greatest possible flow of 
water. Extensive field and laboratory research has been conducted on gravel 
filters, a considerable amount of which applies to wells. The Johnson Well 
Screen Company has done extensive work in the development of natural sedi- 
ments to serve as filters. The size of the slots in the well screen and the 
size and size distribution of the gravel for the artificial envelope should be 
selected principally on the basis of the size characteristics of the sediments 
of the natural formation. 

Many criteria relating to the design of filters have been suggested by 
various investigators. How well these apply to gravel envelopes is not known. 
During recent years, extensive tests to determine the desirable size 
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characteristics for materials for well envelopes have been conducted at 
Colorado A & M College under the direction of Carl Rohwer. Lockman, (7) a 
student of Rohwer’s, reported proposed criteria for filters and well envel- 
opes. A formula proposed by Smith”) is that the ratio of the 50-percent size 
diameters of the pack and the aquifer, (the P - A ratio), be 5. This was con- 
firmed by Leatherwood, (6) who found that instability of the sand occurred for 
values of this ratio exceeding 5.3. These investigators used gravels essen- 
tially uniform in size. Lockman, also using uniform-sized gravels sieved be- 
tween adjacent size sieves, found, nevertheless that the uniformity of the sand 
greatly affected the ability of the pack to resist plugging. He proposed a cri- 
terion involving the uniformity coefficient (ratio of Dg to Dy) of the aquifer 
such that the product of the P - A ratio and the uniformity coefficient be in 
the range 5 to 8. A uniform-sized pack gravel seemed to give best results. 

A comprehensive report of these investigations is presently being prepared 
by Mr. Rohwer and his associates. 


SUMMARY 


The work of Dupuit and the formulas involved from the same general con- 
cept based on assumed steady flow conditions, have provided a valuable basis 
for well engineering and for field testing to determine the permeability of 
underground sediments. Nevertheless, a steady-state flow condition is not 
the usual case, and most well flow comes about by storage depletion. The 
formula proposed by Theis for a non-equilibrium well with subsequent adap- 
tations to special conditions has yielded excellent results and contributed 
greatly to the success of well engineering. Well field developments should be 
designed to deplete storage only to the extent that subsequent periodical re- 
charge may be expected. A development in which production continuously ex- 
ceeds available recharge is doomed to eventual depletion, although in some 
instances the volume of storage is so large that immediate depletion may not 
be imminent. 

A well is, among other things, a mechanism for the transfer of energy to 
an underground body of water. The energy thus transferred causes flow 
toward the well. In inefficient wells, energy losses in and near the well may 
be excessive. Considerable information is available regarding the hydraulics 
of well components, such as screens, gravel envelopes, etc. and of hydro- 
dynamic details existing near the well. This has contributed greatly to the 
art of well construction. Much additional research is needed, however, prin- 
cipally in an effort to better understand hydrodynamic conditions near the 
well and through the well components, and to be able to relate such informa- 
tion to natural geologic conditions and practical well engineering. 
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